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Direct photons provide a insightful tool to study the different stages of a heavy ion 
collision, especially the formation of a quark-gluon plasma, without being influenced 
by the strong interaction and hadronization processes. The yield of direct photons can 
be determined based on the inclusive photon yield and the background from hadronic 
decays. We present a new analysis technique applied to PHENIX Run4 Au+Au dataset. 
It uses strict particle identification(PID) in the Electromagnetic Calorimeter(EMCal) 
and a charged particle veto to extract a clean photon signal. These photons are then 
tagged with EMCal photon candidates with loose PID cuts, which can be reconstructed 
with high efficiency, to determine the fraction of photons originating from ir^ decays. 
Many systematic uncertainties and detector effects cancel in this method. 



1. Introduction 

Recently PHENIX published measurement of direct photons in different centrality 
ranges in ^/snn = 200GeT^ Au+Au colhsions up to pT=14GeV/c the high 
Pt range, direct photons are produced in initial parton-parton scatterings, quark- 
gluon Compton scattering 5 + q — > 7 + g and quark-antiquark annihilation q + q ^ 
7 + 5 which are well described by perturbative Quantum Chromodynamics(pQCD). 
The direct photon yield is determined through the comparison of the inclusive 
photon yield to the expected yield of background photons from hadronic decays. 
The comparison is done in terms of the ratio of N'' jN"^ , or in short 7/7r*' since 
many systematic errors such as energy scale uncertainties cancel. The photon excess 
then can be expressed directly as a double ratio of measured photons per tt" to the 
expectation from hadronic decays, {"i /tt^) measured I {l /T^^^ackground- The double 
ratio in Fig. [1] indicates that there is a significant direct photons signal for pT > 
AGeV/c and it is consistent with the binary scaled pQCD model calculation from 
p+p collisions. Apparently, the direct photons production is in direct contrast to 
high Pt hadrons suppression observed in Au+Au collisions.^ 

At low Pt range where a significant fraction of direct photons is expected to come 
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Fig. 1. Double ratio of mea- 
sured {'y/''r^)Measured invariant yield ratio 
to the background decay {■j/tt'^) Background 
ratio as a function of for minimum bias 
and for five centralities of Au+Au collisions 
at ^/sWn = 200Gey(0-10% is the most cen- 
tral). The solid curves are the ratio of pQCD 
predictions described in the text. 



Fig. 2. Ratio of "/inclusive /ihadronic as a 
function of for minimum bias of Au-l-Au 
collisions at -y/S]vjv = WOGeV. The pur- 
ple triangles are the PHENIX publish result, 
the blue squares are from internal conversion 
method and the black dots are from exter- 
nal conversion method. Color band on each 
data point is the systematic error and the 
extended line is the statistical error. 



from the thermalized medium of deconfined quarks and gluons, the measurement 
is very challenging. These so-called thermal photons carry information about the 
initial temperature of the Quark Gluon Plasma(QGP). But the huge background 
from tt" decay and the relatively large systematic uncertainties at low pT limit our 
ability to make any statement about this direct photon source. 

In order to improve and understand the direct photon production in the low 
and medium px range, PHENIX has already carried out several new analysis tech- 
niques on RUN4 Au+Au ^snn — 200Gey dataset, including using low mass dielec- 
tron pairs from internal conversionJ^I and external photon conversions in the beam 
pipe.^ Unlike conventional direct photon method with the EMCal, these two new 
techniques both take the advantages of the excellent capabilities of the PHENIX 
detector to measure electrons. Due to the excellent resolution of charged particles 
at low momenta, a significant measurement for I < px < 5GeV/c was achieved(see 
Fig. [2]). The method presented here still uses photon candidates from EMCal and 
conventional double ratio technique but instead of measuring tt^ yield in the dou- 
ble ratio of {"//tt") Measured /{i/t^^) Background, wc measurc 7^" which are photons 
from TT^ — > 77 decay by reconstruction of photon pairs invariant mass. In this new 
method we will have to compare 7/7^^ for inclusive photons in data with the cor- 
responding one for hadronic decay photons, which we can obtain from Monte-Carlo 
simulation. Since 7 and 7^^ are basically the same physics quantities, many correc- 
tion factors such as detector efficiency and acceptance cancel out explicitly and so 
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do the systematic uncertainties introduced by these factors. 
2. Analysis 

The data presented in this proceeding were collected during the 2004 Au+Au 
y/IWN = 20006^ run of RHIC by the PHENIX experiment. The EMCal in the two 
central arms {\ri\ < 0.35) is used to obtain the raw inclusive photon candidate. The 
EMCal consists of two subsystems: six sectors of lead-scintillator sandwich calorime- 
ter(PbSc) and two sectors of lead-glass Cherenkov calorimeter(PbGl). In this anal- 
ysis, we only use the PbSc part and events with the vertex \zygrtex\ < SOcm from 
the collision center. That leaves us about 700 million minimum bias events after 
various quality assurances. A series of stringent PID cuts are applied on photon-like 
clusters in the EMCal based on the time-of-flight(TOF) and the shower profile(x^). 
Furthermore, charged particle contamination is removed by associating clusters in 
the EMCal with charged hits in the pad chambers(PC3) positioned directly in front 
of EMCal so it provides a clean photon sample. A minimum px cut of 0.2GeV/c is 
also applied. 

2.1. Decay Photons tagged by tt" 

In each event, these clean photons N-y are combined with photons reconstructed in 
the EMCal under loose PID cuts(x2 & TOF k pr > 0.2GeV/c but without PCS 
veto) and their invariant mass is calculated, generating a foreground of physics 
pairs. The combinatorial background is removed by event-mixing technique, which 
creates uncorrelated photon pairs from different events. The mixed event spec- 
tra is normalized to the corresponding real photon pair invariant mass distribu- 
tions below(50- 100Afe^^/c2) and above(180 - SOOMeV/c^) the tt° peak then sub- 
tracted from the foreground. Fig. [3] shows normalized background together with 
foreground(left) and subtracted foreground(right) for clean photon px range from 
5.0GeV/c to 5.5GeV/c. A Gaussian plus second-order polynomial function is used 
here to fit the 7r° mass peak since there is still some background residual left be- 
neath the peak after subtraction. So the contribution from tt" 77 tagged as N^^o 
can be determined by interpolating the Gaussian fit. 



The number of clean photons and the contribution tagged as coming from 
decay are measured as a function of pr of clean photons. As Eq. [1] indicates, the 
measured yield of clean photons depends on the reconstruction efficiency e and 
the geometry acceptance a. Since charged tracks are removed by the PCS veto cut, a 
hadron contamination correction mostly coming from neutron and antineutron has 
to be applied where Xhadron is the ratio of the remaining hadrons to all clusters 
that satisfy the clean photon's PID cuts. The tagged photon yield A^„o(Eq. [2]) 



N^{Pt) = e X a X N'^^\pt)I{1 - Xuadron) 

No {pt) = e X a X eL X / X (1 - pconv) x N^^lfipr) 



(1) 
(2) 
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Fig. 3. Left panel: Invariant mass of tight-PID-cuts photons and loosc-PID-cuts photons in same 
event and normalized mixed events for tight-PID-cuts photons with 5.0 < Pt < 5.5GeV/c. Black 
points are foreground and red points are normalized combinatorial background. Right panel: In- 
variant mass of tight-PID-cuts photons and loose-PID-cuts photons after background substraction 
for tight-PID-cuts photons with 5.0 < pT < 5.5Gey/c. A fit of Gaussian plus a second order poly- 
nomial is drawn. 

has an additional dependence on the loose-PID-cuts photons' efficiency e^, the 
conversion correction pcov and the conditional probability / to find a loose-PID- 
cuts photon in the PHENIX acceptance, given that its tt" partner clean photon is 
already reconstructed and accepted(see Section [2?2|) . Therefore, by calculating the 
ratio between and N^^o , the efficiency and acceptance for clean photons are 
expHcitly canceled out. 

The loss of N^„o due to photon conversion 7 — s- e+e~ is corrected and the 
conversion correction is estimated by the material budget between the collision 
vertex and the PCS in front of the EMCal. The conversion material in RUN 4 
has been estimated around 5 — 7% /Xq and it is slightly different among different 
detector areas. Since we use a tagging method to measure clean photons from tt" 
decay, only loose-PID-cuts photons are affected due to the conversion loss. Therefore 
a correction factor of (1 — Pconv) ~ 94% is applied. 

2.2. Simulation 

A similar ratio between the hadronic decay photon yield N^hadTo^ (Eq. [3]) and the 
tagged photon yield from decays N^^o (Eq. |4]) can be calculated from simulations 
where the acceptance a is canceled again. 

N^Haaro„ {pt) =ax iV;^Si.„„ (pt) (3) 

iV^,o(pT) = ax/xiV;:S'(pT) (4) 

The decay photon calculations are determined by a fast Monte-Carlo simu- 
lation of 7r° and r/ decays. The PHENIX measured 7r° spectrum^ is used and 
the rj distribution is calculated based on 77^-scaling. The normalization factor 
of ?7/7r" = 0.45 ± 0.10 at high pr is used.lSSlrhe EMCal energy resolution of 
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Se/E = 5% © 9%/VE is applied and a detailed description of the EMCal ac- 
tive areas from data are used in the simulation. From the Eq. (jSj and Eq. Q we 
will have 

{N^H..ro,./N^^o) calculated - (1//) X (iV;^Si™„ /7V™S') (5) 

On the left side of Eq. E (^"hSi™„ /A^'" ) is a straightforward physics quantity 
of excess photons contribution from other hadronic sources than tt" and depends 
only on the rj/n^ ratio and rj — > 77 decay branching ratio. So after applying the 
conditional probability / which is depends on 7r° decay kinematics and the EMCal 
active areas, the ratio {N^hadron /N^„o) calculated can be directly compared to the 
measured one. 

In order to estimate the reconstruction efficiency cl of loose-PID-cuts photons, 
a PISA (PHENIX Integrated Simulation Application) simulation of single particles 
with the complete PHENIX setup based on the GEANT package is conducted. The 
data from the simulation then is merged with the EMCal data from the real events. 
By reconstructing the properties of these embedded particles, the occupancy effect 
of the detector in Au+Au collisions can be estimated. After including all the loose- 
PID-cuts, the efficiency is determined to be around 82% independent of pr beyond 
the minimum pT cutoff of Q.2GeV/c. 

The hadron contamination contribution Xhadron is also calculated with a full 
PISA simulation, tt*^, tt^, and pp are generated as input particles by using 
the actual PHENIX measured spectra.^^lphe spectrum of nn is estimated based on 
measured proton and antiproton cross sections. For clean photons iV^ since most 
of charged tracks are removed by the PCS veto cut, the neutron and antineutron 
contamination is dominant and it is largest around pt of 2GeV/ c due to the contri- 
bution from annihilating antineutrons. Then it is effectively reduced by the shower 
shape cut and becomes negligible above 4G'ey/c, which is less than 1%. 



3. Result 

The comparison of the ratio of measured N^/N^^o and calculated N ^hadron /N^^a 
is essentially the same as the double ratio shown in Fig. [TJ Any significant excess 
of the double ration beyond unity indicates a direct photon signal. 



N^^_ ^ ^ {N-,/N^„o)measured 

TTTnci — X (1 Pconv) X (1 -^hadron) X — — — - 

A^™ad,.o„ [N^h.dro^ /N^^o )calculatcd 



The preliminary RUN4 Au-|-Au results, in Fig. [4l show minimum bias double 
ratio R.y by using this new tt" tagging method compared with RUN2 published 
result and RUN4 preliminary result from conventional method. The main sources 
of systematic errors are the uncertainties in the N^^o yield extraction and the 
hadron contamination correction Xhadron, which gives a final systematic error on 
the double ratio at the level of 10%. All results agree with each other within the error 
bar. For pT > 2GeV/c, the tagging method result offers a significant direct photon 
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Fig. 4. Double ratio R^, as a function of photon pT for minimum bias in Au+Au y^sjviv = SOOGeV 
collisions. The solid black squares are RUN4 tagging method result, the blue triangles are RUN4 
conventional method result(points shifted) and the solid red dots are RUN2 PHENIX published 
result. The color band around data points is the combined error. 



measurement above 1. But for px below 2GeV/c it is still limited by systematic 
errors. Compared with conversion methods mentioned in Section [1] it has smaller 
statistical error and has the ability to extend to mid px range. 



4. Summary 

With the extended statistics provided by the RHIC 2004 Au+Au run and this 
new powerful techniques to extract direct photon in medium and low px range, a 
significant direct photon signal can be measured and results agree with PHENIX 
published result and several other new independent methods including internal 
and external conversion. These promising measurement with improved systematic 
uncertainties indicates that for px < 5GeV/c direct photon yield is lying above 
NLO pQCD expectation and suggests a possible thermal photon emission. 
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